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We present a comprehensive analysis of the performance of different classes of Gaussian states 
in the estimation of Gaussian phase-insensitive dissipative channels. In particular, we investigate 
the optimal estimation of the damping constant and reservoir temperature. We show that, for 
two-mode squeezed vacuum probe states, the quantum-limited accuracy of both parameters can be 
achieved simultaneously. Moreover, we show that for both parameters two-mode squeezed vacuum 
states are more efficient than either coherent, thermal or single-mode squeezed states. This suggests 
that at high energy regimes two-mode squeezed vacuum states are optimal within the Gaussian 
setup. This optimality result indicates a stronger form of compatibility for the estimation of the 
two parameters. Indeed, not only the minimum variance can be achieved at fixed probe states, 
but also the optimal state is common to both parameters. Additionally, we explore numerically 
the performance of non-Gaussian states for particular parameter values to find that maximally 
entangled states within d-dimensional cutoff subspaces [d < 6] perform better than any randomly 
sampled states with similar energy. However, we also find that states with very similar performance 
and energy exist with much less entanglement than the maximally entangled ones. 



PACS numbers: 03.67.Hk, 03.65.Ta, 42.50.Dv 

I. INTRODUCTION 

Decoherence lies at the core of all difficulties in imple- 
menting quantum information technologies. It degrades 
information being transmitted, stored, and processed, in 
an irreversible way. All these processes can be thought 
of different kinds of quantum channels, decoherence in- 
evitably affecting all of them. We arc interested in assess- 
ing the deviation from ideality in Gaussian channels, and 
the precision attainable when tested with Gaussian re- 
sources. Namely, we consider a dissipative thermal bath 
with mean photon number N and damping constant 7, 
and explore how well different Gaussian resources per- 
form in identifying these parameters. Our question is 
immediately relevant to the field of quantum information, 
since assessing the deviation from the identity channel, 
i.e., the ideal information transmitter, is the principal 
requirement to implement large scale quantum commu- 
nication. The need for an efficient characterization of dis- 
sipation in continuous variable systems is becoming a req- 
uisite for a number of quantum information tasks, such 
as quantum repeaters 0, or quantum memories 1 
among others. The burden of dissipation is also hinder- 
ing advances in cavity QED |j| 0] and superconducting 
quantum circuits ||. 

On the other hand, measuring decoherence is not only 
relevant for quantum information technology. In several 
contexts, decoherence can be related to physical quanti- 
ties of practical interest, e.g., photon loss is strongly re- 



lated to impurity doping concentration in semiconductor 
lasers @. Nonlinear magneto-optical effects [l(| can be 
understood as photon loss, an effect with several techno- 
logical applications such as low-field magnetomctry and 
gas density measurements [Tl| - il3j . Additionally, photon- 
photon scattering in vacuum is still an unobserved predic- 
tion of both quantum electrodynamics and non-standard 
models of elementary particles [l4[ . These are only a few 
among the several applications that involve the precise 
determination of losses in dispersive media. For this rea- 
son, we will pose our problem and formulate our results 
in a general theoretic formalism, in order to keep our 
results as general as possible. 

We address the problem of estimating the parame- 
ters of a Gaussian channel describing the dynamics of 
a bosonic mode a, coupled with strength 7 to a thermal 
reservoir with mean photon number N. In the interac- 
tion picture, and within the Markovian approximation 
(at any time the mode and the bath remain unentan- 
gled), the completely positive dynamics and the action 
on the mode a is described by the superoperator 

5(7, AO = cxp 1 (NL[<f] + (N + l)L[a}) , (1) 

with L[o]p = 2opo^ — op — po< o. For convenience we 
arrange the channel parameters in the two-dimensional 
vector = (7, N) and let 9 be the estimator, correspond- 
ing to the outcome of the final measurement. 

Previous related work in the literature addresses the 
problem of estimating a state within a Gaussian family 
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in several different situations. Among others we should 
mention the works of Yuen & Lax [lSj ] . Helstrom [l6j], 
Holevo [l2|, Hayashi fisll . Adesso & Chiribella and 
Hayashi & Matsumoto [20|. However, all of these works 
focus on estimating some of all possible parameters of 
a Gaussian state. Most of them focus on estimating 
either displacement or temperature for fixed degrees of 
squeezing, while others consider the degree of squeezing 
in a vacuum state. No work exists, to the best of our 
knowledge, that addresses the problem of estimating all 
parameters of a Gaussian state. If such a work would 
exist, the problem of estimating a Gaussian channel with 
Gaussian probe states would reduce to a subproblem of 
the former. However, the lack of such a general result 
demands a dedicated solution. 

The setup that we consider is quite generic. We al- 
low for i) extending the channel of interest to include an 
ancillary mode b, unaffected by the channel (i.e., iden- 
tity superoperator T) , obtaining the channel S* = S £g> I, 
ii) choosing any bipartite Gaussian state po, of which k 
copies will be sent through the channel and Hi) perform- 
ing a generalized measurement [characterized by a 
POVM {M {k) }, Ali k) > 0, / d6Mf ] = 1], on the collec- 
tive state = (S*po)® k . This allows for a rather gen- 
eral scheme, which tests the channel with independent 
and identically prepared probes, while the generalized 
collective measurement may include arbitrary quantum 
transformations applied to p g k ^ prior to the measurement 
[sec Fig. [T] . 

As a quantifier of the quality of the estimate one can 
use the covariance matrix 

V e (M) = J d§(9-9)(9-9) T tv[p M § }, (2) 
which may depend on the chosen measurement M = 




FIG. 1: (Color online) Scheme for measuring the values of 
parameters 7 and N. N copies of a bipartite state p are 
independently prepared, and acted upon by k instances of 
the unknown channel I® S(6). A collective measurement A 
is finally performed. 



{Mx} and the particular channel being tested, 9. Al- 
ternatively, as is customary in standard statistical infer- 
ence, we can relate any error cost function 1(9, 9) to the 
covariance matrix by performing a Taylor expansion of 9 
around 9, and defining Gg = \d 2 £(9,9)\g =e as the Hes- 
sian of 1/2, thus 

(£) =tr[G 6 Ve(M)]+o(9-9) 2 . (3) 

More generally, an arbitrary positive semidefinite weight 
matrix Gg can be defined to account for the relevance 
assigned to each parameter. Two extreme cases, where 
one only cares about one or the other parameter can 
be accounted for by the choices G 7 = diag(l,0) and 
Gn = diag(0, 1). This approach allows also to define 
strategies where one is only interested in a particular 
linear combination of the parameters X = x^9^, by 
setting Gx = XX T . With these considerations, and 
given a large number k of copies of the fixed probe state 
po, one can ask what is the smallest possible value of 
(£) = lim /£ ^ 00 fctr[G e V r e (M( fe ))] that is allowed by the 
laws of Quantum Mechanics. 

The main question we wish to answer is the following: 
To what extent can one optimize expression (J3]) by us- 
ing different Gaussian resources? More precisely, we will 
compare the performance of thermal, coherent, single- 
mode squeezed and two-mode squeezed vacuum in esti- 
mating the damping 7 and the temperature N at different 
parameter regimes, comparing the performance of each 
probe state at equal amount of energy input to the chan- 
nel, n = tr[poa^a]- In order to elucidate the role that each 
one of these resources plays in the estimation problem, 
we will make some simplifying assumptions. Namely, we 
will not consider the combination of different resources, 
e.g. two-mode squeezing and displacement. 

A priori it would seem that the problem may have 
several different variants depending on the chosen cost 
function (or G matrix). We will see, however, that some 
general statements can be made. Anticipating the results 
that will be presented in the present work, we will prove 
that: 

1. Choosing a two-mode squeezed vacuum input state 
Po, the parameters 7 and N can be optimally es- 
timated simultaneously. That is, no compromise is 
required in the optimization of Vg(M). This holds 
true even when the optimal measurements for 7 and 
N do not commute. 

2. For both parameters 7 and N, and at any given en- 
ergy, two-mode squeezed states always outperform 
any other class of Gaussian states. 

The combination of these two statements unveals a strong 
compatibility between the problem of estimating the 
damping 7 and the temperature N in the Gaussian set- 
ting. In Section [H] we develop our approach and derive 
the first result of our paper, namely, that optimal pre- 
cision with two-mode squeezed states can be attained 
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simultaneously for both parameters, thus allowing to di- 
vide the problem of computing precision bounds into two 
independent problems. Sections IIIII and IIVI explore the 
precision bounds for estimating 7 and TV individually, 
comparing the performances of different Gaussian re- 
sources, focusing especially on some physically relevant 
regimes of the parameters. In Section |V| we move on 
to explore numerically the non-Gaussian arena and com- 
pare the relative performances with respect to Gaussian 
probes. Section IVT1 concludes the paper with a discussion 
and an overview of the obtained results. Details of the 
technical proofs are provided in two appendices. 



II. THE ULTIMATE QUANTUM LIMITS 

The Heisenberg relations place a fundamental limit on 
the precision with which one can measure any given ob- 
servable. When it comes to quantities not associated to 
an observable, as is our case, it is necessary to resort to 
quantum estimation theory, which studies the fundamen- 
tal quantum mechanical limits to the precision of mea- 
surements in a variety of situations. A first lower bound 
can be obtained from Hclstrom's Fisher information ma- 
trix G1II3, 



V e (M) > J{0)- 



VM 



(4) 



where J (6) = Retr[pg AA T ] is defined as the covariancc 
matrix of the symmetric logarithmic derivatives (SLD) 
A M fulfilling dps/ 89^ = A M op, with AoB = {AB+BA)/2. 
The inherent non-commutativity of Quantum Mechan- 
ics forbids, in general, to attain this inequality when 
the problem is multi-parametric, as in our case. Opti- 
mizing the measurement for one parameter will in gen- 
eral compromise the measurement precision on the oth- 
ers. When considering single-parameter estimation prob- 
lems, it is well known that local adaptive measurements 
attain Eq. ^ [2l], [HI- However, even if the optimal 
measurements for both parameters do not commute, it 
may still be possible to devise a measurement strategy 
to attain simultaneously both bounds. Recent progress 
in the theory of Local Asymptotic Normality for quan- 
tum states [23|, H3| suggests that equality in Eq. Q is 
asymptotically attainable if and only if [251 ] . 



tr[p g [A M , A„] 



0. 



(5) 



The SLD's for our problem, as well as in more gen- 
eral contexts, were obtained by the authors in [261 ] . 
In Appendix [A] we prove Eq. © for the case of two- 
mode squeezed vacuum probe states. The implications 
of Eq. ([5]) are two-fold. On one hand, it allows to prove 
that asymptotically, the estimation problems for both pa- 
rameters become independent, so that they can be ana- 
lyzed separately. This will be the subject of the two next 
Sections. On the other hand, we will show that two- 
mode squeezed states form the optimal Gaussian class of 
states. As a consequence, we will prove the existence of 



and provide the explicit expression of precision bounds 
for the simultaneous estimation of 7 and N for error cost 
functions that have diagonal G matrices. 

The single-parameter precision is quantified by the 
asymptotic standard deviation, 



A6 




6)Hr\ 



(6) 



which is bounded by the quantum Fisher information. 
The latter will be generically denoted as J 7 or Jn de- 
pending whether we are considering the yield for param- 
eter 7 or N , respectively. When no confusion arises, we 
will omit the subscript. When we are referring to a par- 
ticular yield of QFI, specific to a given class of states, we 
will denote it with the corresponding subscript J co h., Jth., 
J sq . and J2_ m for coherent, thermal, single-mode squeezed 
and two-mode squeezed vacuum states, respectively. Pre- 
cision bounds are then given by 



A6 > 



1 



(7) 



We will call Jg the yield, or performance. 

The obtained results should be interpreted in the fol- 
lowing way. The two single-parameter problems have, as 
optimal observables the corresponding SLD's [27j . which 
can be explicitly computed from [26| and are quadratic 
in the creation and annihilation operators. In the mul- 
tiparametric case, with diagonal G matrices, the bounds 
obtained in the following two sections provide all neces- 
sary quantities needed to determine the asymptotic error 
cost. The optimal measurement, however, will require 
a general collective measurement, which is likely to be 
beyond the technical capabilities of present-day technol- 
ogy. We will focus on the theoretical attainable precision 
and not discuss the details of the implementation of the 
optimal observables. This will, nevertheless, provide a 
means to gauge the efficiency of more applied studies 
such as tomographic (28j . single- mode Gaussian [29[ and 
non-Gaussian [30(, or entanglement-assisted schemes 31 



In order to determine the precision attainable with dif- 
ferent Gaussian resources we consider single-mode probes 
parameterized as 



Pa =D(d )S(r )p„ Q S\r )Dt(d ). 



(8) 



where p v oc f ~pjy§ J , is a single-mode thermal state, 
D(d) = cxpi(d 2 Q - d x P) [d = (d 1 ,d 2 )\ and S{r) = 



exp i (ra 



1-2 



"a ). The first moments and covariancc 



matrix of the single-mode states (So) are given by do 
and 



2-<o — — 









3 -2r 



while the energy in the probe is given by 
i^ocosh2r + |d | 2 — 1 



(9) 



(10) 
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Using the results of [26| it is easy to obtain the yield in 
Fisher information as a function of the final parameters 
for both the single- and two-mode probe states. The 
yields for single-mode probes are 



d\e 



-2r 



2 a 2r 



die 



2i> 



+4 [N + 



2 - 1 
; 2 cosh 4r 



^ - 1 



1\ i/cosh2r 



J N = 4(e^ - 1) 



2 1 + v 2 cosh4r 



(11a) 
(lib) 



expressed as functions of the state parameters after the 
action of the channel S. In order to obtain the final values 
for the yield, one needs to consider the three different 
situations, n = |do| 2 /2, 2n+l = vq and 2n+l = cosh2ro, 
and substitute in Eqs. (fTTj) . 

Turning to entangled probe states, we know from [26| 
that the optimal state for constraints of the form 
tr[po o^a] < n can only be pure. Since we are interested 
in evaluating the sensitivity of the different Gaussian re- 
sources, and two-mode squeezing is the only genuinely en- 
tangling resource, we will restrict ourselves to two-mode 
squeezed vacuum states, that we will denote by 



pS = Sam(ro)|0>(0|4 n (r ) J 



(12) 



where |0)(0| denotes the two-mode vacuum. S2m(r) = 
cxp i(ra^ — r*ab) denotes the two- mode squeezing op- 
erator, where b is the ancillary mode. Notice that we 
consistently denote parameters in the probe state with 
the subscript 0. The covariance matrix for the two-mode 
squeezed state is 



/ cosh 2r 



y,* - 



o 



— sinh 2r 







V 



cosh 2r sinh 2r 
■ sinh 2r cosh 2r 
sinh 2r cosh 2r 



(13) 



whereas the energy reads 



1 



-(cosh2r 



1) 



(14) 



for the two-mode entangled state. We will consistently 
use n to compare the performances of the different classes 
of states. 



III. ESTIMATING LOSS 7 

We now proceed to analyze the problem of estimat- 
ing 7 alone, under the assumption that the mean photon 
number N is known. This problem has been partially 
addressed in the literature j29l - l31j with different degrees 
of generality. In previous studies, emphasis is placed 
in zero temperature channels (N — 0). In [3l| several 



distinct probe states arc considered, always with fixed 
tomographic measurements X and P, whereas [2^, [3(| 
focus their attention in the optimal probe states, respec- 
tively within the single-mode Gaussian states and non- 
Gaussian states. Moreover, while 3l| considers the use of 
entangled probes, no consideration is made about the op- 
timality of the measurement scheme. On the other hand, 
References [29|, |30( consider optimality of both measure- 
ment and single-mode probe states, but they do not con- 
sider the use of entangled probes. In this section we will 
combine both approaches, namely, considering different 
Gaussian resources, including entanglement, while still 
using the powerful tools of quantum estimation theory 
in order to take into account the corresponding optimal 
measurement for each probe state. 

The yield for two-mode squeezed probe states as a 
function of the final state parameters is a highly involved 
expression which provides no physical insight. Remark- 
ably, plugging in the dependence of the final parame- 
ters as functions of the initial mean-photon number and 
the channel parameters provides manageable expressions, 
reported in Appendix [B] along with the yields for single- 
mode states. Notice that in both cases, two-mode squeez- 
ing is taken at fixed phase. This does not affect the gen- 
erality of the analysis since single- and two-mode squeez- 
ing along different quadratures can always be taken to 
the standard forms © and (fT3"|) by means of single-mode 
phase shifts. The phase insensitivity of the considered 
channels guarantees that this will not affect the yield in 
Fisher information. 

We thus proceed to compare the different resources by 
writing the output parameters as functions of the initial 
ones [32j , which in turn are functions of the available en- 
ergy. The resulting general expressions are exceedingly 
complex and provide no particular physical insight. Wc 
will, instead, explore specific parameter regimes of phys- 
ical interest. In order to present the forthcoming results 
in a manageable form, we will use, when convenient, the 
following definitions, 



x = n(n + 1), 
y = N(N + 1), 
z = e 7 -l. 



(15a) 
(15b) 
(15c) 



A. Zero-temperature baths, N = 

As a first approach and in order to put our results in 
context with previous related studies [29l - [3ll | we analyze 
the limiting case of baths at zero temperature. Zero- 
temperature baths are the most commonly encountered 
in quantum optics, and to which most of the existing 
literature is dedicated. Under the assumption that N = 
0, the general expressions obtained in Appendix [Bl reduce 
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B. Low energy regime 




FIG. 2: [Color online] Log-Log plot of the yield for the 7 pa- 
rameter using thermal states (red) , coherent states (blue) , sin- 
gle mode squeezed (green) and two-mode squeezed states (ma- 
genta) at different energy regimes, for TV = and 7 = 0.01. 
The saturation of the performance for thermal and squeezed 
states is clear. On the other hand, the linear dependency 
of coherent and two-mode entangled states is readily visible. 
This is in accordance with Eqs. (|16[) . 



to 



Jcoh. — 

Jth. = 

J2-m 



I 

n 



z - 
n 



1 



z 1 + 2(2 + 2(71 + 1))' 
n 



(16a) 
(16b) 

(16c) 
(16d) 



It is easy to verify that J%- 
quantitics. The relations J2-1 



is the largest of all these 
> Jcdh. > Jth.. arc obvious. 



On the other hand the relation J2_ m > J S q. only requires 
to observe that (1 + z 2 )(l + z(z + 2(n + l))) -1 < 1. 

A characteristic feature that we will encounter later on 
with greater generality is the fact that both thermal and 
single-mode squeezed states saturate their performance 
when n is large [See Fig. [2]. On the other hand, the 
performance of coherent and two-mode squeezed states 
grows linearly with n. This clearly leaves the latter as 
the two candidates for optimality when n is large. A 
relevant question is when does the performance of one 
become much larger than the other's? As will be seen, 
two-mode squeezed vacuum states perform always better 
than coherent ones. On the other hand, it is easy to see 
that the increase in performance of two-mode squeezed 
states is most relevant when z<z|l(7Cl), since 
then J2_ m 3> J C oh. ■ This increase in performance is inde- 
pendent of the amount of energy in the probe state. This 
can be observed also from Fig. [2] 



We now turn to the most general situation where the 
thermal bath has nonzero temperature, i.e. photons 
can leak into the quantum system additionally to leak- 
ing out from it. We consider the two interesting pa- 
rameter regimes with practical relevance, namely, that of 
low-energy probes and that of high energy probes. The 
former is best suited for situations where the properties 
of the channel (bath) under inspection are sensitive to 
the effect of intrusive probing. It is worth stressing that 
in some cases, the gain by using a small amount of en- 
ergy may not provide a substantial gain w.r.t. the per- 
formance of the vacuum. On the other hand, there are 
situations where the choice of the probe state critically 
determines the attainable accuracy. We wish to identify 
those situations. 

In the low energy regime (n <C 1) we perform the Tay- 
lor expansion 



J^J^ + J^n + oin 2 ), 



(17) 



where obviously is independent of the kind of state 
being considered, and corresponds to the performance of 
the vacuum. The common leading constant is thus given 
by 



J(0) 



N/i 



l + z(JV + l) 
Corrections of order n contribute with coefficients 



J 



(i) 



1 



coh. 



y th. 



J (1) 



J. 



(1) 



2-m 



l + z(l + 2/V)' 
(z + l)(l + 2z(/V+l)) 
z 2 (l + z(A + l)) 2 ' 
2N+1 z(l + N) 2 (l + 2N) 
I (l + z(/V + l)) 2 

2(l + 2/V) 2 
+ (1 + z) 2 + 2Nz{\ + z(N + 1)) 
(z + 1) 2 + N{z(z + 2) + 2) 
z(l + z{N+l)) 2 



(18) 

(19a) 
(19b) 

(19c) 
(19d) 



where we have defined y = N(N + 1). A few comments 
are in order. First of all, notice from Eq. 



(|B8a)l that the 
yield for coherent states is a polynomial of first degree 
m n. Therefore, + j££ n gives the exact expression. 
On the other hand, the thermal correction J±P < is 
always negative, which implies that weak thermal fields 
perform worse than the vacuum [See Fig. [3] (Left)]. 

The first question to ask is when do coherent states 
provide any significant improvement over the vacuum. 
For this we consider the condition j' ) <C J^L n. This 
reduces to 



N{z(2N + !) + !) 
z(z(N + + 



< n. 



(20) 
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FIG. 3: [Color online] Yield for the 7 parameter using thermal states (red), coherent states (blue), single mode squeezed 
(green) and two-mode squeezed states (magenta) at different energy regimes, for N — 0.9 and 7 = 0.3. [Right] Log-log plot 
of the yield for n values greater than 1. The linear behavior of coherent and two-mode squeezed states is readily apparent, 
whereas saturation of the yield occurs both for thermal and single-mode squeezed states. At this parameter values we have 
(e 7 — 1)(2N + 1) ~ 0.98, thus J2- m /J co h. — 2.0 for n S> 1. [Left] Detail of the yield at low n values in linear scale. Different 
slopes corresponding to the different values of j' 1 ' are apparent. 



On the other hand, the same condition for two-mode 



squeezed probe states <C J^n, reduces to 



N(z(N +!) + !) 
(N + l){z + l) 2 + 2N 



< n. 



(21) 



Notice that for moderate values of N and large values of 
7 (z > 1) Eqs. ([201) and (HJ) reduce to N(2N+l)/(z(N+ 
l))<n and N/z <C n respectively, which shows that the 
regime required for Eq. (J2TJ) is entered earlier than that 
of Eq. (|2H| for increasing z. 

In the limit of small losses {z <C 1), we can expand 
Eqs. (|2TH) and (f2~Tj) to zcroth order in z to obtain, for 



N ,r2 

— + N 2 < 71 

z 



(22) 



whereas for <^ J^m n we S e t 



N 



1 + 2N 



< n. 



(23) 



The condition for -C J^-m 71 thus reduces to N <C 
nj (1 — 2n) ~ n. Notice that this is not a sufficient condi- 
tion to achieve a significant improvement using coherent 
probes [Eq. ([22])]. 

On the other hand, one can see that 7 <§C 1 (that 
is, z < 1) with moderate TV (z(l + 2iV) < 1) implies 

"^coh ^ This means that at moderate tempera- 

tures, and with low energy in the probe state, two-mode 
squeezed states significantly outperform coherent states 
in the regime of small losses. 



C. High energy regime 

High energy probes are of interest when the channel 
being probed is not as delicate, or needs not be pre- 
served. A natural instance of this situation is in prob- 
ing the photon-photon scattering predicted by QED and 
non-standard models of elementary particles. The high 
energy regime has a substantially different behavior. Ex- 
panding the relevant yield functions from Appendix [Bl in 
inverse powers of n we obtain series of the form 



J 7 = J<- 1) n + j(°)+o(l/n) 



(24) 



In some cases J^ -1 -* will vanish, rendering the corre- 
sponding class of states useless compared to those for 
which J^ 1 ) does not vanish. Explicitly, we have 



J coh. 
■7th. 

J2-m 



l + 2 (2iV+l) 



o(i), 



1 + [ - 

n 



z(2N + 1) 



0(1) 



(25a) 

(25b) 

(25c) 
(25d) 



The first relevant fact to notice is that, contrary to the 
low energy regime, different Gaussian resources perform 
differently in the limit of large energy This is hardly a 
surprise. We observe that the asymptotic performance 
is bounded for thermal and single-mode squeezed states. 
This is in contrast to the fact that single-mode squeezed 
states have proven highly efficient for other precision 
measurements such as optical phase {33U35| and mag- 
netometry [3fj| . among others. On the other hand, co- 
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FIG. 4: [Color online] Yield for the N parameter using thermal states (red), coherent states (blue), single mode squeezed 
(green) and two-mode squeezed states (magenta) at different energy regimes, for N = 0.9 and 7 = 0.3. [Above] Detail of the 
yield at low n values. The X parameter is negative, which implies that single-mode squeezing is detrimental (at small n values). 
[Below] Behavior at high n values. 



herent and two-mode squeezed states provide an un- 
bounded yield, as is manifest by the nonvanishing j'" 1 ' 
terms, giving a linear growth in J 7 with increasing n. 
This makes the constant zeroth order correction irrel- 
evant. A very important difference between coherent 
states and two-mode squeezed states becomes readily ap- 
parent. While for coherent states, the rate of growth of 
is bounded, for two-mode squeezed states it is not. 
In particular, the difference between the two yields be- 
comes most significant when z(2N + 1) <C 1 and negligi- 
ble when z(2N + 1) 3> 1. This is certainly relevant for 
detecting very small damping parameters, for which the 
inverse dependence in z may even be sufficient to over- 
come the practical limitations of achieving very high n 
values. 



the coupling constant cannot be chosen at will. On the 
other hand, this does not necessarily provide the optimal 
estimation accuracy. As we will see, two-mode squeezed 
states can outperform the sensitivity of the vacuum state 
or other classes of Gaussian states, for any value of the 
parameters. 

Contrary to the effect of the decay parameter 7, that 
affects both first and second moments, the temperature 
in the bath only affects the second moments. First mo- 
ments evolve independently of the bath temperature [32| . 
This has immediate consequences for the sensitivity of 
coherent states, which will perform equivalently to the 
vacuum state. Thus, we will not consider coherent probe 
states in this section. We will follow the same approach 
taken in the previous section, by addressing different en- 
ergy regimes in the probe states. 



IV. ESTIMATING TEMPERATURE 



Quantum thermometry has become a subject of high 
physical relevance with the advent of ultracold atomic 
gases [37M39I ] . At low temperatures, new methods need 
to be envisaged to determine the magnitude of thermal 
fluctuations in atomic clouds. In this section we analyze 
the quantum-limited precision bounds to the estimation 
of temperature (mean photon number N) in a bosonic 
thermal bath, coupled to a probe system prepared in a 
Gaussian state, with the coupling strength not necessar- 
ily large, i.e., far from thermalization. 

A straightforward method to measure temperature is 
to let the probe system coupled to the bath to thermalize. 
Quadrature measurements then provide an estimator of 
the mean photon number, providing in turn an estima- 
tor of the bath temperature. This approach has several 
drawbacks. Most importantly, it requires in general a 
large coupling constant 7 » 1, in order to reach the 
steady state. However, there may be situations where 



A. Low energy regime 



As in the case for J 7 , Jn will be limited to a constant 
value corresponding to the sensitivity of the vacuum state 
J*- ' at small values of n. However, since first moments 
are unaffected by the temperature of the bath, the yield 
of coherent states will equal that of the vacuum. As in 
the previous section, we will expand Jn in powers of n, 
in order to obtain and analyze the low-energy yield for 
each class of states, 



J. 



N 



J(°) + J«n + o(n 2 ) 



Taking the expressions from Appendix [B] we get 

z(z + l) 2 



N(l + z{N + l)) 



(26) 



(27) 
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FIG. 5: [Color online] Yield for the N parameter using thermal states (red), coherent states (blue), single mode squeezed 
(green) and two-mode squeezed states (magenta) at different energy regimes, for N = 0.7 and 7 = 0.08. [Left] Detail of the 
yield at low n values. The X parameter is positive, which can be see by the positive slope of the single-mode squeezing (n 1). 
[Right] Behavior at high n values. 



to which corrections of order n contribute with factors 
(z + l) 2 (z(2/V+l) + l) 



J, 



(i) 



Hi 



J. 



J, 



N 2 (z(N + 1) + 1) 2 ' (28a) 
(1) = 32z(z + l) 2 (2(| - 1) - z (4zg 3 + (z + 2)| + l)) 
sq - (4ze + 4£ - z - 2) 2 (z (4z£ 2 + 4£ + z + 2) + 2) 

(28b) 

(!) (2/V+l)z(z + l) 2 



2-m 



N(N + l)(z(N + l) + iy 



(28c) 



where we have defined £ = N + 1/2. 

It is immediate to observe that < which im 
plies that, similarly to the situation for J 7 , small ther- 
mal fluctuations in the probe state can only be detri- 
mental. On the other hand, the small n correction for 
single-mode squeezed states has no unambiguously de- 
fined sign, the latter being positive only when X = 
2(£ - 1) - z (4z£ 3 + (z + 2)£ + 1) > 0. Solving the in- 
equality X > for z we obtain 



z < 



(2g-l)y^+l-2g-l 
+ 1) 



(29) 



The right hand side is positive only when £ > 1, which 
means that only for N > 1/2 it is possible to have a 
gain in yield by using single-mode squeezed probes. We 
thus conclude that for < N < 1/2 the best single- 
mode probe (at low energies, n < 1) is the vacuum 
state. However, if the bath temperature is sufficiently 
high [TV > 1/2], it is possible to improve the sensitivity 
of the vacuum by single-mode squeezing. 

Wc now turn to analyze the yield of two-mode squeezed 
states. As can be readily seen in Figs. [4] and \5\ the behav- 
ior is rather simple. The slope at small values of n is al- 
ways positive and, moreover, as follows from the relations 
in Appendix [Bj always larger than that of single-mode 
squeezed states. Imposing that J%~.n » yields 



greater than unity, we have n(27V+l) 3> 1 and, given that 
n is small, it follows that N has to be large. Therefore 
we can reduce the condition to z < (2n — 1)/N, which 
is never satisfied for small n. We thus conclude that 
low energy entangled probes always outperform single- 
mode ones, for any regime of the parameters, although 
the improvement cannot always be of significant mag- 
nitude (J2-m n cann °t be achieved for small n 
values). Most remarkably, two-mode squeezed states are 
the only class of Gaussian states that perform better than 
the vacuum for any values of the parameters. 



B. High energy regime 

The high energy limit for Jjy is quite uninteresting. 
The reason being that all yields saturate, and no sig- 
nificant improvement can be achieved using two-mode 
squeezed probes as compared to the use of coherent 
states. The limiting expressions read, to order 1/n 



Ji 



th. 



0. 



J: 



2-m 



Z 2 (l +zf 

n 2 
2(1 + z) 2 
(2iV+l) 2 (2N + l) 3 zn 
(1 + z) 2 {l + z) 2 



(l + zf 



N(N+1) N (2N 2 + 37V + 1) zn 



(30) 
(31) 
(32) 



We have presented the second order term for J t h. be- 
cause it is the leading one. The yield, however, tends 
to vanish, as can be expected by observing that in the 
limit of a highly energetic probe, the thermal fluctuations 
in the probe are infinitely larger than those induced by 
the bath, and therefore no inference about the latter can 
be obtained. Single- and two-mode squeezed states have 
nonzero limiting yields, but, as anticipated, the yield sat- 
urates for highly energetic probes. It is easy to check 



n(2JV+l) ^> z(N+l) +N+1. Since the right hand side is that in the limit n — > 00 we have J2- 



> J sq . but no 



9 



order can be established between the vacuum state and 
single-mode squeezed states. A crossover can occur be- 
tween the performances for the vacuum and single-mode 
squeezed states depending on the values of 7 and N. 

Concentrating on the asymptotic yield, a natural ques- 
tion is to understand in what situations does two-mode 
squeezing perform much better than the vacuum state. 
This is simple to answer by observing that 



lim 



J, 



1 



1 



coh. 



(N + l)z 



(33) 



so that if (N + l)z <C 1 (in the high n limit) then 
^2-m 3> J C oh.- This condition is likely to occur in several 
practical situations, for channels close to ideal, i.e., when 
the decay rate 7 is very low and the temperature is low or 
moderate. Notice that this is not the first time that we 
encounter this condition for a significant improvement 
of squeezed states over coherent ones. All results indi- 
cate that whenever there is an improvement of two-mode 
squeezed states over the vacuum, it is always in the limit 
of small coupling to the bath (7 <C 1, i.e., a close to ideal 
channel) . 



V. ENTANGLEMENT AND THE 
PERFORMANCE OF NON-GAUSSIAN STATES 

So far we have seen that, in order to optimally estimate 
the channel parameters in Eq. ((T|), two- mode squeezing 
is the most effective resource within the arena of Gaus- 
sian inputs. This fact suggests two immediate questions. 
1) Is there a direct relation between entanglement and 
the performance in channel estimation? And 2) Are 
there non-Gaussian states outperforming the squeezed 
vacuum for the same energy supply? Questions simi- 
lar to 1 ) have arisen previously in the literature, both in 
the finite dimensional case ptJ - l42j | and in the continuous 
variable one (albeit in somewhat different setups [43]]). 
considering different kinds of channels. Question 2) has 
been also addressed in the context of quantum chan- 
nel estimation [30l |. and regarding the performance in 
other quantum information tasks, especially continuous- 
variable quantum teleportation [44| . 

With the available techniques it is difficult to give a 
precise quantitative answer, since numerical techniques 
are not well-suited to deal with infinite dimensional sys- 
tems, and analytic control methods are not yet well de- 
veloped beyond the Gaussian regime. The difficulty with 
numerical methods resides in the fact that Hilbert-space 
truncation is rendered useless at N 7^ 0, because ther- 
mal baths immediately populate all levels in Fock space, 
thus rendering a direct numerical approach futile, or, in 
the best of all cases, a rude approximation. In order to 
address the question of sensitivity of non-Gaussian states 
for parameter estimation, techniques beyond those devel- 
oped so far are needed, and it is beyond the scope of this 
work to pursue them. Instead, we will provide a simple 
qualitative answer by restricting the channels of interest 



to those at zero temperature, i.e.. N ~ 0. In this case, 
populated levels in Fock space only decay to lower en- 
ergy levels and Hilbert-space truncation provides exact 
numerical results. 

We have computed numerically the quantum Fisher 
information J 7 at 7 = 0.1 for 4000 states picked from 
C 4 (£> C 4 constituting the subspace with at most 3 pho- 
tons in each mode, randomly distributed according to the 
SU (16) Haar measure. In Fig. [5] we report our findings. 
The left plot displays the values of J 7 against the mean 
photon number in mode a, and the right plot displays 
the efficiency (i.e., the ratio J 7 /n) against the entropy 
of entanglement E(p) = — trp a logp a , p a = titPo- Along 
with the random states we display the results for maxi- 
mally entangled states within the cutoff 3 < d < 6 and 
the two-mode squeezed vacuum. Our findings reveal that 
most states have a very high performance in relation to 
the amount of energy they have in the channel mode. In 
particular, for a fixed dimension cutoff the more entan- 
gled the states are, the better they perform on average. 
In particular, the maximally entangled state attains the 
bound set by the squeezed vacuum. It is of course inter- 
esting to ask how does the performance change when one 
introduces temperature in the channel. This is beyond 
the capabilities of numerical methods relying on dimen- 
sional cutoff, and more advanced methods would need to 
be envisaged. 



VI. DISCUSSION 

We have obtained the sharp precision bounds on es- 
timation of 7 and N for four classes of Gaussian states, 
namely, coherent, thermal, single-mode squeezed vacuum 
and two-mode squeezed vacuum. We have shown that 
the two-mode squeezed vacuum always outperforms any 
other class of Gaussian states. The improvement of two- 
mode squeezed vacuum states versus coherent states is 
most relevant when the coupling parameter 7 is weak. 
In particular, at zero or finite temperature, the yield J 7 
of two-mode squeezed states increases much faster with n 
than the yield of coherent states at small values of 7. For 
Jjv, comparing the yield of two-mode squeezed vacuum 
states versus the yield of the vacuum, we find that, for 
small values of n, no significant improvement can be ob- 
tained. The situation changes dramatically in the high 
energy regime, where, despite a saturation of the yield 
(all yields saturate to a maximum value), saturation oc- 
curs at much higher yields when z(N + 1) <C 1. Summa- 
rizing, we have shown that two-mode squeezed vacuum 
always outperforms any other Gaussian resource for both 
estimation problems (7 and N), and we have identified 
the situations in which this improvement is most signifi- 
cant. We have also provided numerical evidence that the 
squeezed vacuum state provides an upper bound to the 
performance of arbitrary states with dimensional cutoff. 
This suggests a deeper analysis, and the causes for this 
optimality should be investigated. 
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FIG. 6: [Color online] Scatter plot of 4000 random probe states with at most 3 photons in each mode (dots), maximally 
entangled states at dimension cutoff 3 < d < 6, \tp) = Ylk ( crosses )) an d two-mode squeezed vacuum states (solid 

line). Left: Yield against mean photon number in the channel mode a as a reference. The dashed line interpolates the behavior 
of maximally entangled states. Color code corresponds to the entropy of entanglement in the probe state E(p) — —trp a logp a , 
Pa = titp. Right: Ratio between the yield and the mean photon number against the entropy of entanglement. Observe that 
several highly efficient probes (those with high ratio J 7 /n) are relatively unentangled as compared to the maximally entangled 
ones with similar performance. Clearly, the squeezed vacuum state is much more entangled than randomly sampled states with 
similar yield. 



Turning to the multiparametric problem of simultane- 
ously estimating both parameters of the channel, we have 
shown that the optimal state (two-mode squeezed vac- 
uum) provides optimal sensitivity for both parameters, 
and that in the many-copy limit (k — > oo) a collective 
measurement exists which saturates both bounds simul- 
taneously, as a consequence of Eq. ([5]). These results 
imply that the problems of estimating 7 and N enjoy a 
strong form of compatibility. Not only the optimal probe 
state is common to both problems, but also the corre- 
sponding optimal measurements commute in the asymp- 
totic limit. Therefore an overall protocol optimizing both 
tasks can be envisaged (in the asymptotic k — > 00 limit). 
It is a relevant question to determine whether lifting the 
restriction to Gaussian states can provide an increased 
performance, and whether such kind of compatibility re- 
mains true in the more general setting involving arbitrary 
probe states. 
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channel p = <S*po- R l = (Qi, Pi, Q2, P2) are the canon- 
ical operators with [R 1 ,R 1 ] = ifl^ and f2 is the 2-mode 
symplectic matrix, f2 = uj © w, with 



1 

-1 



(Al) 



It is convenient to work with the centered canonical oper- 
ators R 1 = R l — (R 1 ) with which we define the covariancc 
matrix of the output state E 1 - 7 = (R 1 o RP). Here we have 
introduced the symmetric product AoB = (AB+BA)/2. 

We will be using Einstein's summation convention and 

'i 



tensors of the form A\ B 1 l will be written as [A 



B] lJ 



y - 



v <g> 



Analogously A 1 \v 3 will be written [A ® v 
if no confusion arises. Finally, we construct a covari- 
ance matrix with covariant-contravariant transformation 
rules, namely E l fc = (R 1 o Rj)fljk so that E = Ef2. 

The covariancc matrix of the input state will be de- 
noted Eq, 



So =2 



/ coshr sinhr 

cosh r — sinh r 

sinhr coshr 

\ — sinh r cosh r 



(A2) 



and observe that the final covariance matrix is in stan- 
dard form 49!, 



Appendix A: Proof of Eq. (O for two-mode squeezed 
vacuum states 



al cZ 
cZ bt 



(A3) 



In this section we show that, for two- mode squeezed 
vacuum probe states, Eq. ([5]) holds. Since we are prob- 
ing channel S* with bipartite states we will be working 
with two bosonic modes. Start by introducing some no- 
tation. Let (X) = tr[pX] be the expectation value of 
any operator, evaluated with the state output from the 



where Z = diag(l, —1), a = e~ 7 coshr + (1 — e~ 7 )(iV + 
1/2), b = coshr and c = a -7 / 2 sinhr. Next, define the 
superoperator 2? M as V^S* = d^S* and notice that, for 
both parameters 7 and N it holds that [2(| 



?V = ",,.- « (/'"/'/'" - {R 3 R l ) ° p) 



(A4) 
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where the a matrices have the diagonal block structure 

(A5) 



a. 



fN+1/2 -i/2 
i/2 TV +1/2 







V 



a N = (e 7 - 1) 







1 



1 







(A6) 







With this, the SLD's read [2J 

A M = afl ,ijC ij (A7) 

where 

= [^^"i + [£ (J) ] i i | fl*ofl' 1 (A8) 

and 

L (0) = ifl-^EgS- j(E®l-l(8)S)), (A9) 
L (1) = 2i(XT 1 ® (i?) - (i?) (A10) 
L (2) = iz) _1 (l<g)l-i(E(8)l-l(8)S)), (All) 



and wc have defined 



D = E® E - -1® 1. 
4 



(A12) 



Also, observe that both a M matrices share the same struc- 
ture = k^P + u^Q where 



P = 



o o J ' w loo 



(A13) 



In order to evaluate Eq. (J5]) , observe that some terms in 
Eq. (|A8p immediately drop when taking the commutator 
or the expectation value 



tr[p[A M , A„]] = a^a^i'j'T^^ , 



(A14) 



where 

rpiji'j' 



tv[p[C ij ,C ij }} (A15) 

= [ i(i) ] i i[ i(i) ] i '^ fefe ' ( Aie ) 

+ [i( 2 )]^[L( 2 )]%([i? fc o R l ,R k ' oR 1 '}). 

Proceed by evaluating T yi . The following relations 
will be used, 

tr[p[R l ,R 1 '}} = ifi«', (A17) 
\x\p[&,i#' oF]] =0, (A18) 

tr^ir o ff' o =i(n" , E 3 ' :7 " + q^'e^" 

+ ^"E J ' 4 ' +0«"E"). (A19) 



Let us define the quartic momenta K 1 ^ 



i j 



R l ,R k o B 1 ]}, and the superoperator (in phase-space) 
D as D[A]ij = AitjiD % \-. Notice that D is linear and 
25 [A T ] = T) [A] T , thus T>[A] preserves the symmetric or 
antisymmetric character of A, and so does any function 
of ID, such as ID -1 . With this we de fine V = D~ 1 [P] and 
Q = D~ 1 [Q], and wc can write Eq. (|A14j) as 



,, ,, .,.,t % ^i 



kk 



(A20a) 



(A20b) 

+ i{n^ v - ^)PijQi>j'KiAfk,f>K klk ' 1 ' (A20c) 



Concentrating on line (|A20a[) . it is easy to see that it 
trivially vanishes. Due to the antisymmetry of we 
have 



'J 

k 



(^Pij+i^Qi^L^Yl (A21) 
= i^QijlL^Yl. (A22) 

Then (|A20aP reads 

(Qv[LVf k ) (QML^,)n kk ' =0(A23) 



due to antisymmetry of il. Line (|A20b[) also vanishes 
identically due to the antisymmetry of K and the sym- 
metry of {KftKvVijVv j' - i^iyQij Qi> j'Wki-^k'i' u nder in - 
terchange of (kl) O (k'V). We are left with line ([AlOc)) . 
Note that 

rklk'V 



Va Q. vy MslM/J'vK klk ' 1 ' (A24) 



where we have replaced M by their symmetric [Ms — 
1 ® 1] and antisymmetric components [Ma = *(E <£> 1 — 
1 (g) E)] in accordance with the symmetry of V and Q 
respectively. Using Eq. (|A19|) we can rewrite Eq. (|A24[> 



as 



VijQi'j'M^M(,{,K klk ' 1 ' = 4 (te[QEETO] + tr[QE7>E] 

(A25) 

By virtue of tr[AB] = tr[A T B T ], Q T = -Q and 
(E"PE) T = E"PE the second term trivially vanishes. We 
are thus left with 

tr[p[A M , A„]] = 4i(* M i„ - ^^)tr[QEOEPf2]. (A26) 

In order to prove that this vanishes, let us analyze the 
structure of Q and V. Let us express £> -1 L4] as 



2-^ = -4 



^4 fc (E T ) fc AE* 



.fc=0 



(A27) 
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Given any matrix A in block form 
1 c Z 60 1 



(A28) 



we can define A k = (E T ) fc AE fc , and it is easy to show 
that 



E T A fc E = 



a k+1 l Cfc+iZ 
c k+ iZ b k+ il 



(A29) 



with cifc+i = a fe a 2 + c(cb k + 2ac k ), b k+ i = a k c 2 + b(bb k + 



2cc k ) and c fc+ i 



aafcC — bb k c — abc k . From 



this it is straightforward to solve the recurrence equations 
and perform the resummation in Eq (|A27[I to obtain the 
explicit form of V. Similarly, one can obtain the explicit 
form for Q. For our purposes, however, it is only relevant 
to notice that 



v _ 1 a p t c p Z 
1 c p Z bp! 



a q io c q t 

— C q l bqLU 



(A30) 
(A31) 



-c p (a q + C q )Z Cp(c q - b q )t 
On the other hand we have 



. (a 2 - c 2 )uj -c(a - b)X 
1 c(a-6)X {b 2 -c 2 )uo 



Finally, the trace in Eq. (|A26[1 reads 



Finally, we can obtain the structure of VflQ, which reads 
-(a p a q + c p c q )t (cLpC q - c p b q )Z \ (A32) 



1 

1 

(A33) 



tr[QEQETO] = -(a p a q + c p c q )(a 2 - c 2 )tr[Z] 
+(a p c q — c p b q )c(a — 6)tr[ZX] 
+c p (a q + c q )c(a — 6)tr[ZX] 
+c p (c q - b q )(b 2 - c 2 )tr[w] 
= 0. (A34) 



This shows that 



taMA^AJ] =0. 



(A35) 



To recapitulate, the expectation value of the commuta- 
tor of the SLD operators vanishes for probe states given 
by the covariance matrix of Eq. (|A2p . In order to sec 
that this result applies to all possible two-mode squeezed 
vacuum states, one has to observe that if Eq = SJ^qS t , 
where S are local phase shifts. Hence, the phase invari- 
ance of the channel guarantees that the output covariance 
matrix transforms in the same way, E' = SE5 T , thus one 
can write 

tr[p[A M , K]] = a tl , ij a v , kl S i i ,S j j ,S k k ,S l l ,T^' k ' 1 ', (A36) 

where the left hand side refers to the state corresponding 
to E' and T refers to the state E. Finally, since S is a 



symplectic orthogonal transformation, and a M is a linear 
combination of the symplectic and the Euclidean metrics, 
these are invariant under the congruence, 



ap ltij S % i ,S : ' j , = a 



fti, i'j' ■ 



(A37) 



This completes the proof of Eq. ©, i.e., that the 
Heisenbcrg limit for 7 and N corresponding to two-mode 
SVS can be, asymptotically, simultaneously achieved. 



Appendix B: Exact expressions 

We report, for completeness, the exact expressions for 
the Fisher information for both parameters 7 and N. 
When convenient, we will use the variables 



x = n(n + 1), 
y = N(N+1), 
z = e? - 1. 



(Bl) 
(B2) 
(B3) 



Moreover, we find it useful to define t = n + N + 2nN = 
x + y — A 2 where A 2 = (n — N) 2 . In terms of x and y 
we have 



VI + 4x - 1 



N 



(B4) 
(B5) 



t = vd+ ^d+^-i (B6) 



The following inequalities are trivial consequences of the 
previous definitions, 



x > 0, 

y> 0, 

z > 0, 

x + y > t> 0. 



(B7a) 
(B7b) 
(B7c) 
(B7d) 



Working out the results of [26[ we obtain, for the dif- 
ferent classes of states, 



X 



N/z 



7 " l + z(N + l) l + z(2A+l)' 



(B8a) 



J^ q - = — f 2 ^ + V) v (B 8b) 

7 z t + yz (z + l) 2 + 2z{t + yzY K ' 

jth. = x + y-t 



J, 



yz 2 + tz + x 

2-m _ t + XZ 



z((t+i)z + iy 



Sc) 
(B8d) 
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For N we get 

^■coh. 



z(z + l) 2 
(l + z(N + l))N' 

TSq. 

N 



(B9a) 



J 



th. 



7 2-m 

N 



[z(z + 1) 2 (1 + 4x + z(2i + 2yz + z + 2))] 

x [i + 2(ar(8y + 2) + y(z(2yz + 2 + 2) + 3)) 

+te 2 (l + 4y)], (B9b) 

z 2 (z + l) 2 , . 

(B9c) 

(B9d) 



x + z(t + yz) ' 
(t+l)z(z+l) 2 



+ z + 1) 
It can be checked that 



T2-m 
J 7 


> 


rth. 


(BlOa) 


r2-m 
J 7 


> 


TCOll. 

1 


(BlOb) 


T-2-m 
J 7 


> 


TSq. 
J 7 ■ 


(BlOc) 



Although apparently nontrivial, a systematic method can 
be used to check all the above inequalities. We illustrate 



the procedure by considering Eq. (IBlOap . We first write 
it as 



t + xz 



x + y — t 
z((t + l)z + 1) ~ yz 2 + tz + x 



> 



(Bll) 



to obtain 



(t + xz) (yz 2 + tz + x) - z(x + y- t)((t + l)z + 1) > 0. 

(B12) 

Now arranging in powers of z we obtain 

xyz 3 + (t(t + 1) - x - y) z 2 

+ (x 2 +t(t + 1) - x - y) z + tx > 0. (B13) 

Now it is trivial to show that all coefficients are positive 
by observing that t(t + 1) = x + y + 4xy. The same 
method applies to all inequalities. Notice that since 7 
always appears as exp 7, then only integer powers of z will 
appear in the expressions. No series expansions will be 
necessary and therefore one only needs to check positivity 
for a finite number of coefficients. 
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